Kinetic and thermodynamic study on the esterification of oleic acid over SO3H-functionalized eucalyptus tree bark biochar catalyst

Herein, esterification of oleic acid (OA) over tosylic acid functionalized eucalyptus bark biochar (TsOH-MBC) to synthesize fatty acid methyl ester (FAME) was investigated. The TsOH-MBC catalyst was prepared via pyrolysis-activation-sulfonation process at various impregnation ratios and was characterized by SEM, FTIR, EDX, XRD, BET, TGA and acid site density techniques. The catalytic performance of the sulfonated biochar catalyst was described in terms of acidity and FAME yield. 6 g of sulfonic acid loaded on 10 g of MBC (6TsOH-MBC) appeared to be most appropriate combination to achieve a highly active catalyst for the esterification of OA with 96.28% conversion to FAME at 80 °C for 5 h with catalyst loading of 4.0 wt% and 8:1 methanol/OA molar ratio. The catalytic reaction kinetic data were very well described by the second-order model, with a rate coefficient of 0.223 mL mol−1 h−1 at 80 °C and activation energy of 81.77 kJ mol−1. The thermodynamic parameters such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta H$$\end{document}ΔH, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta S$$\end{document}ΔS and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta G$$\end{document}ΔG were determined to be 78.94 kJ mol−1, 135.3 J mol−1 K−1 and 33.03 kJ mol−1, respectively. This research provided an environmentally friendly procedure for FAME production that could be replicated on a commercial scale.

activated charcoal to biochar, however, it is necessary to improve the latter's physicochemical properties. The chemical treatment of biochar with zinc chloride (ZnCl 2 ) as a modifying agent can result in activated biochar with a hydrophobic surface, a high surface area, a well-developed porous structure, and an abundance of surface functional groups 1,13 . Yusuff et al. 14 discovered that activating poultry litter biochar with ZnCl 2 produced an adsorbent with a high sorption capacity for crystal violet removal from aqueous solution. Thus, ZnCl 2 -activated biochar was considered as a structure support for the development of a heterogeneous acid catalyst.
It has been reported that the synthesis of SO 3 H (sulfonic acid group) functionalized carbon based material via sulfonation using H 2 SO 4 produces a catalyst with high acidic strength and better activity 1,15,16 . Dechakhumwat et al. 3 reported that a solid acid catalyst prepared by impregnating corncob biochar with H 2 SO 4 increased the yield of esterified product. However, using H 2 SO 4 as a sulfonating agent for sulfonated catalyst development can cause corrosion if a high pressure withstanding reactor is not used, necessitating special precautions 11,17 . Thus, there is a need to find alternative sulfonating agents that are less corrosive, have a low vapour pressure and have better catalytic activity. As a result, tosylic acid (TsOH) or p-Toluenesulfonic acid (p-TSA) has been discovered to be a suitable replacement for H 2 SO 4 because the former possesses the aforementioned properties 3,18 . Russo et al. 19 used a SO 3 H-functionalized C/SiO 2 composite obtained through TsOH and H 2 SO 4 sulfonation for biomass conversion and discovered that TsOH prevented damaging corrosion.
Several studies on the kinetics of fatty acid esterification over various catalysts have been conducted. Many of the reported studies advocated for the use of a pseudo-first-order kinetic model to evaluate the kinetic parameters [20][21][22][23][24] . A kinetic study with a second-order model for esterification of oleic acid over zinc acetate catalyst had been reported, and it agreed reasonably well with the experimental data obtained 25 . Nonetheless, there is a scarcity of information on the kinetics and thermodynamic analysis of oleic acid esterification over a biomass-derived solid acid catalyst.
The current research focused on the synthesis of a SO 3 H-functionalized eucalyptus tree bark biochar catalyst for oleic acid esterification. To investigate the physicochemical properties of the as-prepared catalyst, several characterization techniques (SEM, FTIR, EDX, XRD, TGA, BET and acidity) were used. Kinetics and thermodynamic analyses of the catalytic reaction process were performed using different model equations. In addition, the stability test was carried out under optimal reaction conditions to determine the possibility of catalyst reuse.
Catalyst preparation. Preparation of ZnCl 2 modified biochar. Firstly, the collected ETB was thoroughly washed with clean water to remove sand, dried at 80 °C for 6 h in an oven, and then ground into powder using mortar and pestle. After that, ETB powder was sieved through 0.3 mm sieve mesh to achieve the needed particle size. To make biochar, the sieved ETB powder was pyrolyzed in a muffle furnace with no oxygen (at 500 °C for 1.5 h and a ramping rate of 5 °C/min). The pyrolyzed material (biochar, BC) was allowed to cool before being stored in a glass container. To activate the biochar, 10.0 g of BC powder was mixed with 200 mL of 3 M ZnCl 2 solution in a beaker, stirred on a magnetic stirrer (for 5 h at 60 °C) and filtered through filter paper. The activated biochar was then washed several times with distilled water to remove impurities until the pH of the solution reached neutral. After that, the sample was dried at 110 °C overnight. The obtained chemically activated biochars are henceforth referred to as ZnCl 2 -modified biochar (MBC).

Preparation of SO 3 H-functionalized eucalyptus tree bark biochar catalyst.
A fixed amount of the MBC sample (10 g) was suspended in distilled water, and different amounts of TsOH (2, 4, 6, and 8 g) were gently added and stirred for 5 h at 60 °C on a magnetic stirrer. The resulting homogenized solution was then oven dried for 24 h at 100 °C. The dried sulfonated biochar catalysts were denoted as mTsOH-MBC, where m denotes the mass of TsOH used for MBC sulfonation, e.g., 2TsOH-MBC denoted that 2 g of TsOH was used for the sulfonation process.
Analysis of MBC and TsOH-MBC catalyst. The surface functional groups, phase composition and crystallographic structure, thermal decomposition trend, surface morphology and elemental composition and textural characteristics of the MBC and optimal mTsOH-MBC catalyst were examined by using Fourier transform infrared (FTIR) spectrometer (Perkin-Elmer Spectrum TWO, spectra range 4000-400 cm −1 ), powder X-ray diffractometer (XRD) (D8 Advance, Bruker AXS Gmbh, Germany), TGA/DTG analyzer (Shimadzu DTG60, Japan; analysis condition: N 2 carrier flow rate = 100 mL min −1 , heating rate = 15 °C min −1 and range of temperature = 30-900 °C), field-emission scanning electron microscope (FE-SEM, Quanta 200F, Eindhoven, Netherland), Micromeritics model analyzer (ASAP 2010, USA; analysis condition: the sample was degassed at 250 °C for 5 h to remove adsorbed molecules from the surface, and Brunauer-Emmett-Teller (BET) surface measurement was done under N 2 gas at 77 K), respectively. Moreover, acid density (or acidity) of the optimal sulfonated biochar catalyst was determined using ion-exchange titration procedures reported by Nda-Umar et al. 11 . For the determination of the property, 0.05 g of the optimal mTsOH-MBC catalyst sample was suspended in 20 mL of 2 M NaCl solution and mixed for 1 h on a magnetic stirrer. The suspension was thereafter filtered, and the filtrate obtained was titrated against 0.05 M NaOH using phenolphthalein as an indicator. The acidity, expressed as mmol of H + /g cat, was estimated using Eq. (1). Catalytic activity studies. Esterification of oleic acid. Esterification of oleic acid (OA) with methanol over TsOH-MBC catalyst was conducted in a 100 mL two-neck round bottom flask connected with a reflux condenser and a magnetic stirrer. Throughout the catalytic reaction process, the flask was inserted in a silicon oil bath to control the temperature. For each experiment, a certain mass of the catalyst sample was added to methanol, mixed for 10 min, and 10 g of OA was added to the mixture. After that, the reaction mixture was heated to the desired reaction temperature and continuous stirring began immediately at 500 rpm in order to avoid mass transfer limitation. The reaction process was carried out at 60 °C for 2 h with methanol/OA molar ratio of 8:1 using 2 wt% catalyst loading. Following the completion of the reaction, the reaction products were centrifuged at 7500 rpm for 10 min to facilitate the removal of spent catalyst. To remove unreacted methanol and water, the centrifuged liquid was evaporated using a rotary evaporator at 105 °C.
Esterified product analysis. The methyl ester content in the evaporated product was analyzed by gas chromatographic device (Agilent GC, 7890A, USA) coupled with flame ionization detector (FID) (GC-FID) and a capillary column (J & W DH-5HT; dimension: 15 m × 0.32 mm × 0.1 μL). Helium and methyl heptadecanoate were used as carrier gas and internal standard, respectively. The operating condition of the GC-FID was maintained at 350 °C and 30 °C min −1 for 10 min. The FAME yield (Y) and conversion to FAME (X) were calculated using Eqs. (2) and (3), respectively.
where A indicates overall peak areas of FAME from C 14 to C 22 , A I is the peak area associated with the internal standard, V I and C I are the volume and concentration of internal standard, respectively, m is the sample weight, M F is the mass of FAME obtained and M OA is the mass of oleic acid used.
Kinetics and thermodynamic studies. The kinetics and thermodynamic behaviours of OA esterification over 6TsOH-MBC catalyst was studied at different values of temperature (55, 65 and 80 °C) and time (1-5 h). The kinetics and thermodynamic parameters were evaluated as explained below: Investigation of the reaction parameters. The strategy used in estimating the esterification reaction parameters (rate constant and reaction order) has been previously reported by Toor et al. 26 . Esterification of OA with methanol, a reversible catalytic reaction process, is illustrated stoichiometrically as follows: Because the amount of methanol (M) used exceeded the stoichiometric amount, the irreversible non-elementary reaction rate equation is as follows: Expressing Eq. (5) in terms of conversion, we have Inserting Eqs. (6a-c) in Eq. (5) and rearranging, we have: where C OA o , C M o t, k , α , β , C OA , C M , X OA , ∅ are initial concentration of oleic acid (mol mL −1 ), initial concentration of methanol (mol mL −1 ), reaction time (h), rate constant (unit depends on reaction order), reaction order of oleic acid, reaction order of methanol, final concentration of oleic acid, final concentration of methanol, fractional conversion of oleic acid and methanol/oleic acid molar ratio. www.nature.com/scientificreports/ Since many studies on esterification of organic acid with alcohol reported either zero or first order 24,25 , it was presumed in this study that the order of reaction between OA and methanol could be pseudo-first or secondorder, thus indicating that the overall reaction order must be equal to or less than two ( α + β ≤ 2).
Case 1: Assuming pseudo-first-order reaction rate equation by considering α = 1 and β = 0, thus Eq. (7) upon integration becomes: where k 1 is the rate constant for first-order reaction (h −1 ) Case 2: A second-order reaction rate was considered by assuming that α = 1 and β = 1. Therefore, Eq. (7) upon resolution by partial fraction, integration and rearrangement becomes: where k 2 is the rate constant for second-order reaction (mL mol −1 h −1 ).
Estimation of thermodynamic properties. By fitting the experimental data at different temperatures to the Arrhenius equation (Eq. 10), the activation energy (E a ) and pre-exponential factor (A) can be calculated.
A plot of Ink against 1 T yields a slope = − E a R and intercept = InA. Furthermore, because the esterification reaction process was carried out at different temperatures, the thermodynamic properties of OA esterification over 6TsOH-MBC catalyst, including enthalpy change ( H ), Gibbs free energy (( G ) and entropy change ( S ), were estimated using the following expressions [27][28][29] .
Approval and compliance with regulation. This research work was formally authorized by the Academic and Research Unit, Afe Babalola University, Ado-Ekiti, Nigeria. Experimental research and field studies on collected eucalyptus tree bark comply with relevant institutional guidelines and legislation.
Results and discussion FAME yields of the catalyst samples. In this study, four samples of TsOH sulfonated MBC were synthesized, corresponding to the four different masses of TsOH used during the sulfonation process. TsOH-MBC catalysts were obtained in the following concentrations: 2TsOH-MBC, 4TsOH-MBC, 6TsOH-MBC and 8TsOH-MBC. To assess the activity of each catalyst sample, they were used to catalyze the esterification of oleic acid with methanol, and the values of FAME yields obtained were plotted against the different TsOH-MBC catalyst samples used, as shown in Fig. 1.
As shown in the results (Fig. 1), the MBC appeared to be ineffective for converting OA to FAME due to a lack of acidic sites on its surface, preventing the reaction from reaching equilibrium. However, when TsOH was loaded onto MBC, there was a significant increase in FAME yield with increasing the amount of the sulfonating agent used in the sulfonation process. The trend shown in Fig. 1 indicated an increase in FAME content from 32.6 to 92.3%, with a corresponding increase in TsOH mass from 2 to 6 g. When 8 g of TsOH was loaded on the same amount of MBC, the FAME yield decreased significantly. The possible reason for the decreasing of catalyst activity with increasing TsOH amount is the TsOH penetration through the MBC in various ways depending on the uptake level of catalyst support towards sulfonating agent. The TsOH was well dispersed on the MBC surface at low concentration of TsOH, but at higher concentration of TsOH, the acid molecules deposited on the surface of MBC and blocked the active sites 1,30 . These findings indicated that 6TsOH-MBC was the most effective of the prepared catalysts, so it was chosen for further investigation.
The sulfonation process, in general, introduces both sulfonic acid (-SO 3 H) and oxygenated acid (-COOH) sites. The -SO 3 H functional groups, according to Zhang et al. 31 , act as the primary catalytic sites, while the -COOH groups promote the inherent activity of the -SO 3 H, facilitating oleic acid esterification. Protonation of methanol molecules occurs infrequently during the esterification reaction due to the strong nature of the -SO 3 H SEM analysis. Figure 2 depicts the morphologies of the MBC and 6TsOH-MBC catalyst samples. MBC had a sponge-like, rough and irregular surface structure, as shown in Fig. 2a, with many pores on its surface. However, as a result of the impregnation effect, the pores were no longer visible after TsOH was loaded on MBC (see Fig. 2b), indicating even dispersion of the SO 3 H group on the support. The FTIR results (see Fig. 3) corroborated this, revealing a transformation in the chemical structure of the catalyst support, with the formation of oxygen-containing functional groups such as SO 3 H, C-O-H, and OH groups observed. Furthermore, the SEM image revealed that the crystallinity of the 6TsOH-MBC was more pronounced than that of MBC, which was supported by the XRD pattern (Fig. 4), which revealed that the former had a more obvious amorphous carbon structure. It is worth noting that the cracks observed on the solid acid catalyst were most likely caused by partial collapse of the porous structure during the sulfonation process, as reported by Ngaosuwan et al. 1 . Nonetheless, the reacting molecules were able to easily penetrate through the available pores on the catalyst surface, as      www.nature.com/scientificreports/ compared, it was possible to conclude that the sulfonic acid group was successfully inserted into the biochar, implying that acidic sites were active centers during the reaction process.
XRD analysis. X-ray diffractograms of the MBC and 6TsOH-MBC catalyst samples are shown in Fig. 4, with the former revealing some peaks that were attributed to ZnO phase at 2 θ around 34°, 36° and 56° probably owing to the ZnCl 2 used as a modifying agent. Ngaosuwan et al. 1 reported similar finding whereby ZnO phase was detected in XRD diffractogram of ZnCl 2 -activated coffee residue biochar. However, a visible change in the XRD pattern was noticed after sulfonation of MBC sample, as seen in the pattern of 6TsOH-MBC. A peak at 2 θ = 27.1° was noticed in the diffractogram of the solid acid catalyst, which corresponded to the (002) plane of amorphous carbon. This result was in line with the findings reported by Zhao et al. 36 and Dechakhumwat et al. 3 who synthesized sulfonated biochar catalysts from waste pomelo peel and corncob, respectively. Table 1, the MBC was made up of mainly C followed by O with traces of Zn and Cl owing to activating agent (ZnCl 2 ) used for activated biochar, whereas 6TsOH-MBC catalyst contained C, O, Si, S, Cl and Zn. The presence of sulfur in the catalyst sample suggested successful impregnation of -SO 3 H group on the modified biochar. The composition of C and O in the 6TsOH-MBC catalyst decreased probably due to the sulfonation effect. In terms of C/O ratio, the sulfonated eucalyptus biochar catalyst exhibited larger ratio as compared to MBC which indicated that impregnation of TsOH on MBC increased oxygen-containing functional groups and contributed to the retention of oxygen content 13 . More so, it enhanced the hydrophilicity of the as-synthesized catalyst, thereby improving the reactant adsorption onto its surface 1,36 .

EDX analysis. As listed in
Thermal stability analysis. In a bid to investigate the thermal decomposition trends of the activated biochar and sulfonated biochar catalyst, TGA analysis was conducted on both samples. As evident in Fig. 5, the support and solid acid catalyst exhibited different decomposition trends owing to sulfonation effect. Two weight loss stages were observed in the TGA curve for MBC. However, 6TsOH-MBC catalyst exhibited three stages of mass loss, as seen in Fig. 5, indicating that the latter had better thermal stability than the former. The first stage of decomposition of the heterogeneous acid catalyst occurred between 30 and 405 °C, which was as a result of the simultaneous removal of water molecules and decomposition of sulfonic acid (-SO 3 H groups) 11,36,37 . Since degradation of lignin usually occurs at 200-500 °C 38 , the second degradation stage, which took place from 405 to 606 °C,  www.nature.com/scientificreports/ suggested that sulfonation led to complete decomposition of lignin as well as complete removal of hemicellulose and cellulose 39 . However, the last stage of weight loss by the sulfonated biochar catalyst occurred above 600 °C which suggested that the 6TsOH-MBC was thermally stable due to carbon moiety 11,40 , as reaffirmed by the FTIR and EDX analyses.
Textural characteristics and the acidity. The N 2 adsorption isotherms of the MBC and 6TsOH-MBC samples are shown in Fig. 6, with the former exhibiting a type II isotherm and the latter exhibiting a type IV isotherm (mesoporous material). Table 2 presented the textural characteristics of the MBC and 6TsOH-MBC samples. The BET specific surface area, total pore volume and average pore size diameter reduced after sulfonation process. Similarly, Zhang et al. 41 reported a significant decrease in surface area during TsOH impregnation on glucose. The possibility of -SO 3 H multilayer formation on the MBC surface during sulfonation could not be ruled out. The partial blockage of the pores by multilayer dispersion of -SO 3 H reduced the specific surface area 30 , which corresponded to the SEM images of the sulfonated biochar catalyst (see Fig. 2). However, considering the result presented in Table 2, the low surface area of the 6TsOH-MBC catalyst did not signify low activity. The good activity exhibited by the sulfonated catalyst during esterification reaction was certainly due to the presence of acidic sites, considering that the acid density of the MBC was relatively low, 0.02 ± 0.053 mmol of H + /g cat. According to Refaat 42 , catalytic performance is closely related to the acidic sites density and hence the solid catalyst with higher acid density is expected to show better activity during reaction. Thus, this suggested the reason why 6TsOH-MBC catalyst exhibited good performance in esterification of OA.

Esterification reaction condition.
Influence of catalyst dosage. The results obtained during the investigation of catalyst loading influence on the esterification of OA at 65 °C for 5 h with an 8:1 methanol/OA molar ratio using a 1-5 wt% 6TsOH-MBC catalyst dosage are shown in Fig. 7. The FAME content increased with the amount of 6TsOH-MBC catalyst used up to 4.0 wt% and then decreased as the catalyst amount increased to 5.0 wt% due to the increasing viscosity of the reaction mixture, which inhibited the reaction 30,43 . This finding suggested that increasing the catalyst loading could not guarantee increase in FAME yield but raises the overall production cost of biodiesel. As a result, the optimum catalyst dosage of 4.0 wt% was chosen and used for the subsequent esterification reaction studies.
Influence of reaction temperature and time. The influence of reaction temperature on the esterification process at different reaction durations (1-5 h) is illustrated in Fig. 8. The reaction tmperature was investigated at 55, 65 and 80 °C. The esterification reaction conducted at 80 °C for 5 h gave the maximum FAME content of 96.28%,  Reusability of 6TsOH-MBC catalyst. After each experiment, the spent catalyst was collected, washed with n-hexane to remove adhered fatty acid molecules, dried overnight in an oven at 100 °C, and then reused for another esterification reaction. FAME yields of 80.6%, 62.1%, and 38.7% were obtained for three consecutive cycles of oleic acid esterification at 80 °C for 5 h with a methanol/OA molar ratio of 8:1 and a 4.0 wt% regenerated 6TsOH-MBC dosage. This observation was attributed to a decline in the number of acidic sites owing to poisoning of the catalyst surface by fatty acid molecules and washing during catalyst regeneration process, which reduced the FAME content 1,43 . This is supported by the FTIR analysis of the regenerated catalyst ( Fig. 9), which revealed a peak associated with C=O (ester) stretching from fatty acid methyl ester 45,46 . Furthermore, the intensity of peaks (1200 cm −1 and 1040 cm −1 ) associated with the sulfonic group (-SO 3 H) decreased, confirming the loss of acidic sites. This result suggested leaching of the active ingredient (SO 3 H), so a leaching analysis was performed by suspending the solid acid catalyst in methanol and heating the suspension at 55 °C for 2 h. Following that, the methanol was recovered by filtration and used for esterification of OA at 80 °C for 5 h with a methanol/OA molar ratio of 8:1 and a catalyst loading of 4.0 wt%. When the used methanol reacted with OA, a FAME yield of 10.07% was obtained, indicating that leaching of the 6TsOH-MBC occurred. As a result, the decrease in FAME content observed during several runs of the reaction process was caused by catalyst leaching.   Esterification kinetics and thermodynamic studies. Kinetics studies. Figure 10a, which illustrates the linear plot of -In(1 − X OA ) against reaction time at different reaction temperatures (PFO), facilitates the estimation of reaction rate constant ( k 1 ) from the slope of the plot and coefficient of determination ( R 2 ). A linear plot with high R 2 value indicated that such a kinetic model equation is suitable in predicting the catalytic reaction process 30,47 . In this study, however, PFO did not adequately predict the esterification of OA owing to low R 2 values (0.8465-0.9334) for all reaction temperatures, as seen in Table 3.
In an attempt to establish the appropriate kinetic equation for predicting the esterification of OA over 6TsOH-MBC catalyst, the second-order equation was employed and the constants ( k 2 and R 2 ) were estimated at different temperatures from the plot shown in Fig. 10b. As contained in Table 3, the R 2 values were determined to be in the range 0.8798-0.9895 and were greater than those of the PFO equation. Thus, the esterification reaction using 6TsOH-MBC catalyst appeared to be second-order with respect to oleic acid and methanol. Moreover, as also seen in Table 3, the reaction rate constants increased with increasing reaction temperature. The presumed reason was that when the reaction temperature was increased, the collision of molecules was enhanced, which led to decline in mass transfer limitation thus rate of reaction increased 30,48 .
Thermodynamic studies. Because the esterification reaction process was studied at various temperatures, estimating activation energy was critical in studying the effect of temperature on reaction rate 49,50 . The slope and intercept of the plot of Ink against 1/T (Fig. 11) were used to estimate the values of E a and A, which were then used to calculate the reaction thermodynamic parameters such as H , G and S . Their estimated values are presented in Table 4. The values of E a and A were found to be 81.77 kJ mol −1 and 2.45 × 10 11 h −1 , respectively. The activation energy for the esterification of OA over 6TsOH-MBC increased the reaction rate, which was consistent with the value range (33.6-84 kJ mol −1 ) for esterification of fatty acid using solid acid catalyst 51 , indicating that the catalytic reaction was much easier to carry out using SO 3 H-functionalized biochar catalyst as compared to other solid acid catalysts.
The values of thermodynamic parameters for esterification of OA with methanol over 6TsOH-MBC catalyst are also presented in Table 4. The positive values of H suggested that the catalytic reaction process was endothermic. The endothermic and nonspontaneous esterification process was confirmed because the G value for each temperature was positive 48,52 . Moreover, the increase in values of G as the temperature increased suggested an enhanced reaction rate at higher temperature 28 . Entropy, which measures the degree of disorderliness in a reaction system, was also calculated. Table 4 shows that the values of entropy change (�S) for the three reaction temperatures investigated were all positive, indicating an increase in irregularity and randomness at the solid-liquid interface during the esterification reaction. Based on these findings, it is possible to conclude that the 6TsOH-MBC catalyst demonstrated a faster reaction rate, moderate activation energy and mass transfer resistance during the esterification reaction process.

Conclusions
This study found that the SO 3 H-loaded eucalyptus biochar catalyst was effective for the esterification of OA to produce methyl ester. The well-developed porous structure, dominance of active acidic sites and well dispersed TsOH, as validated by SEM, FTIR, acidity analysis and EDX, contributed to the 6TsOH-MBC catalyst performance in the esterification process. At 80 °C for 5 h with an 8:1 methanol/OA molar ratio and 4.0 wt% catalyst loading, a maximum FAME content of 96.28% was obtained. Analysis of kinetic data using two different kinetic models showed that the esterification of OA over 6TsOH-MBC was second-order with respect to each reactant. Studies on thermodynamic behaviour of the reaction indicated that the esterification of OA with methanol over 6TsOH-MBC catalyst was nonspontaneous and endothermic. These findings confirmed that the SO 3 H functionalized eucalyptus biochar catalyst could be widely used for large-scale sustainable biodiesel production.